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I .  INTRODUCTION 


Several  recent  reports1  4  have  described  progress  made  over  the  past 
several  years  in  an  effort  to  develop  a  fully  two-dimensional,  two-phase  flow 
interior  ballistic  model.  This  work  has  focussed  on  the  process  of 
f  amespread  through  the  propelling  charge  as  a  hydrodynamic  problem  and  on  the 
influence  of  the  path  of  flamespread  on  the  formation  of  potentially  dangerous 
pressure  waves  in  the  gun  chamber.  The  effort  was,  to  a  large  extent, 
motivated  by  the  fact  that  early  successes  in  simulating  these  phenomena  in 
Navy  cased-ammunition  guns  '  with  one-dimensional,  two-phase  flow  models  were 
not  reproduced  when  Army  bagged-charge  artillery  became  the  subject  of  study* 2 3 * 5 * 7. 
The  presence  of  circumferential  ullage  external  to  the  bag  apparently  offered, 
at  least  during  the  very  early  stages  of  the  interior  ballistic  cycle,  a 
region  of  high  permeability  capable  of  altering  the  flame  path  and 
equilibrating  longitudinal  pressure  gradients  —  a  process  totally  outside  the 
scope  of  the  one-dimensional  representation.  A  subsequent  quasi-two- 
dimensional  treatment8  recognizing  this  possibility  rapidly  led  the  way  to  a 

fully  two-dimensional  representation  known  as  TDNOVA,  the  subject  of  this 
report . 


’.S.  Gough,  "Two-Dimensional  Convective  Flamespreading  in  Packed  Beds  of 
Granular  Propellant ,  ARBRL-CR-00404 ,  USA  ARRADCOM ,  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD,  July  1979.  (ADHA075326) 

2 

A.W.  Horst  and  P.S.  Gough,  "Modeling  Ignition  and  Flamespread  Phenomena  in 
Bagged  Artillery  Charges,"  ARBRL-TR-02263 ,  USA  ARRADCOM ,  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  MD,  September  1980.  (ADttA091790) 

3 

D.S.  Gough,  "A  Two-Dimensional  Model  of  the  Interior  Ballistics  of  Bagged 
Artillery  Charges,  ARBRL-CR-00452 ,  USA  ARRADCOM ,  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  MD,  April  1981. 

\.U.  Horst,  F.U.  Robbins,  and  P.S.  Gough,  "A  Two-Dimensional ,  Two-Phase  Flow 
Simulation  of  Ignition,  Flamespread,  and  Pressure-Wave  Phenomena  in  the  155- 
MM  Howitzer,  ARBRL-TR-  ,  USA  ARRADCOM,  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD,  (not  yet  released). 

5A.U.  Horst,  T.C.  Smith,  and  S.E.  Mitchell ,  "Key  Design  Parameters  in 
Controlling  Gun-Environment  Pressure-Wave  Phenomena  -  Theory  Versus 
Experiment ,  13th  JANNAF  Combustion  Meeting ,  CPIA  Publication  273,  Vol.  1, 
pp.  341-368,  December  1975. 


A.W.  Horst  and  P.S.  Gough,  "Influence  of  Propellant  Packaging  on 
of  Navy  Case  Gun  Ammunition,"  Jsmr.ma.1  o£  &£L±±iS.ti£.S.  -  Vol.  1 ,  Ho 
258,  1977. 


Performance 
3,  pp.  229- 


A.U.  Horst,  C.U.  Nelson,  and  I.W.  May,  "Flame  Spreading  in  Granular 
Propellant  Beds:  A  Diagnostic  Comparison  of  Theory  to  Experiment,"  AIAA  Paper 
N°-  77-856,  AIAA/SAE  13th  Propulsion  Conference ,  July  1977. 

Q 

P.S.  Gough,  "Theoretical  Study  of  Two-Phase  Flow  Associated  with  Granular  Bag 
Charges,  ARBRL-CR-0038 1 ,  USA  ARRADCOM ,  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  MD,  September  1978.  (ADMA062144) 
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II.  TECHNICAL  DISCUSSION 


a.  Description  of  TDNQVA 

The  TDNOVA  code  provides  an  unsteady,  two-dimensional,  axisymmetric,  two- 
phase  flow  representation  of  the  interior  ballistic  cycle.  As  mentioned 
previously,  the  development  of  TDNOVA  was  undertaken  largely  in  response  to 
the  conf igural  complexities  associated  with  the  use  of  bagged  artillery 
charges,  such  as  the  155-mm,  M203  Propelling  Charge  depicted  in  Figure  1. 
Flamespread  through  bagged  charges  is  believed  to  be  strongly  influenced  by 
the  details  of  the  ullage  which  initially  surrounds  the  bag  and  by  the 
behavior  of  the  bag  material  itself.  Accordingly,  an  explicit  representation 
is  made  in  TDNOVA  of  the  two-phase  region  occupied  by  the  propelling  charge  at 
any  time.  The  flow  in  the  ullage,  which  surrounds  the  region  occupied  by  the 
propellant,  is  treated  as  unsteady,  inviscid,  and  single  phase. 


4,0  OZ.  CLM$  1 
SLACK!  POWDER 
{GEARHART  0WEN5, 


RESIN  IMPREGNATED, 
CLOTH,  VISCOSE-RAYON. 
CLASS  2  (6.5  OZ./YD1) 


JACKET,  LACING 
ACRYLIC-VISCOSE 
RAYON.  CLASS  3 
(DYED  SCARLET,  L„ 

PROPELL  ANT  M3B A  T 


BASE  IGNITER  PAD 
1  OZ.  CLASS  1 
BLACK  POWDER 
(GEARHART  OWENS) 


DONUT 

FLASH  REDUCER 
16  OZ. 
POTASSIUM 
SULFATE 


TUBE  IGNITER 

(MOLDED  NITROCELLULOSE) 
1.4  MAX  O.D.  X  27”  LONG 


LEAD  FOIL  5.5  OZ. 
(DECOPPERING  AGENT) 


(WEAR  REDUCING  ADDITIVE) 
Tl  0j/WAX*1  7.5  OZ. 


Figure  1.  155-mm,  M203  Propelling  Charge 


The  ullage  is  divided  into  several  disjoint  regions,  coupled  to  one 
another  and  to  the  two-phase  flow  in  the  propelling  charge  by  means  of  finite 
jump  conditions  at  all  their  mutual  boundaries.  By  formulating  the  theory  in 
such  a  manner  as  to  use  directly  the  jump  conditions  at  the  boundary,  a 
mechanism  is  provided  for  representation  of  the  influence  of  the  bag. 
Impermeability  is  reflected  directly  within  the  momentum  jump  condition  as  a 
quasi-steady  flow  loss.  Similarly,  the  influence  of  exothermically  or 
endothermically  reactive  components,  such  as  igniter  basepads,  centercore 
tubes,  or  wear-reducing  liners,  may  be  reflected  by  means  of  source  terms  in 
the  mass  and  energy  jump  conditions. 
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The  division  of  ullage  into  the  several  regions  depicted  in  Figure  2  is 
based  on  the  instantaneous  configurations  of  the  external  boundaries  (e.g., 
spindle  face,  chamber/tube  sidewall,  and  projectile  base)  and  the  propelling 
charge  itself.  Each  region  of  ullage  may  be  treated  as  lumped  parameter, 
quasi-one-dimensional  (i.e.,  one-dimensional-with-area-change ) ,  or  two- 
dimensional,  in  accordance  with  user-definable  criteria  based  on  physical 
dimensions. 

As  suggested  earlier,  representation  of  a  basepad  igniter  or  centercore 
tube  may  be  treated  within  the  structure  of  the  bag.  A  centercore  ignition 
charge,  coaxial  with  the  bag,  may  also  be  included  in  the  representation  as  a 
quasi-one-dimensional,  two-phase  flow,  coupled  to  the  state  of  the  flow  within 
the  bag  and  any  ullage  present  at  the  ends  of  the  chamber  by  reference  to 
finite  jump  conditions.  Representation  of  the  ignition  train  also  admits 
specification  of  an  externally  injected  stimulus  of  predetermined  flow  rate 
and  energy. 

Initially,  TDNOVA  provides  a  fully  two-dimensional  analysis  of  flow 
within  the  two-phase  region  occupied  by  the  propellant  bed.  However,  in  all 
calculations  performed  to  date,  the  regions  of  ullage  contiguous  to  the  bag 
have  been  treated  as  quasi-one-dimensional,  the  continuum  coordinate  being 
defined  by  the  common  boundary.  Corner  regions  of  ullage  are  then  given  a 
lumped-parameter  treatment.  Figure  3  illustrates  this  level  of 
representation . 

Following  the  completion  of  flamespread,  rupture  of  the  bag  sidewall,  and 
equilibration  of  the  radial  structure  of  the  pressure  field  to  within  some 
user-specified  limit,  a  quasi-two-dimensional  approach  is  introduced,  similar 
to  that  reported  previously8.  For  the  duration  of  the  ballistic  cycle,  the 
propelling  charge  is  given  a  quasi-one-dimensional  representation,  as  is  the 
circumferential  ullage,  while  regions  of  axial  ullage  at  the  ends  of  the 
chamber  are  treated  as  lumped  parameter  (see  Figure  4). 

Each  of  these  regions  of  continuous  flow  is  mapped  onto  a  regular  figure, 
a  unit  line  or  square,  by  means  of  a  boundary-fitted-mesh-transformation 
algorithm.  The  method  of  solution  is  then  based  on  an  explicit,  two-step 
marching  scheme  which  utilizes  characteristic  forms  of  the  balance  equations 
at  both  external  and  internal  boundaries.  A  detailed  description  of  the  code 
has  been  provided  by  Gough8.  The  reader  is  further  directed  to  an  earlier 
discussion  of  application  of  TDNOVA  to  the  155-mm,  M203  Propellinq  Charqe  by 
Horst  et  al4. 


In  the  sections  that  follow,  we  provide  a  description  of  several  series 
of  baseline  calculations  performed  to  assist  in  evaluation  of  the  operational 
capabilities  of  TDNOVA. 

B.  Camper  1  son  Hith  NOVA 


A  direct  comparison  was  made  between  results  predicted  by  TDNOVA  and  by 
its  quasi-one-dimensional  predecessor  NOVA^.  An  appropriate  data  base,  a 
dimensionally  degenerate  representation  of  the  previously  described  bagged- 


P — — - 

P.S.  Gough ,  "The  NOVA  Code:  A  User's  Manual .  Volume  I.  Description  and  Use," 
IHCR  80-8,  Naval  Ordnance  Station ,  Indian  Head,  ML,  30  December  1980 . 
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Figure  2.  Computational  Regions  of  TDNOVA 


Figure  3.  Two-Dimensional  Mesh  of  TDNOVA 
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charge  problem,  was  devised  by  suppressing  the  presence  of  both  centercore 
igniter  and  circumferential  ullage.  Bag  sidewall  and  centercore  tube 
characteristics  thus  played  no  role  in  the  problem.  Further,  the  basepad  was 
replaced  by  a  predetermined,  one-dimensional  ignition  stimulus.  The  resulting 
problem  is  depicted  schematically  in  Figure  5.  A  summary  of  required  input 
data  for  both  NOVA  and  TDNOVA  is  provided  in  Appendix  A. 

It  is  noted  that,  in  accordance  with  criteria  based  on  size  and  structure 
of  the  flow,  the  NOVA  code  may  assign  a  continuum  representation,  in  the  axial 
direction,  to  regions  of  ullage  at  either  end  of  the  chamber.  Calculations 
were  performed  using  the  NOVA  code  both  unaltered  and  with  a  modification 
introduced  to  maintain  a  lumped-parameter  representation  for  these  regions, 
similar  to  that  provided  by  TDNOVA  after  transformation  to  the  quasi-two- 
dimensional  representation.  In  both  cases,  the  NOVA  simulations  employed  30 
axial  stations,  while  a  30x7  mesh  was  used  for  the  TDNOVA  runs.  Finally,  a 
TDNOVA  calculation  was  performed  with  the  quasi-two-dimensional  representation 
invoked  from  time  zero.  In  the  absence  of  circumferential  ullage,  this  last 
TDNOVA  treatment  was  geometrically  equivalent  to  that  of  the  modified  NOVA 
code . 


A  summary  of  results  is  provided  in  Table  1.  All  values  for  maximum 
chamber  pressure  fall  within  2%  of  one  another,  while  those  for  muzzle 
velocity  fall  within  1%.  Values  for  -AP^,  the  initial  reverse  pressure 
difference  between  breech  and  forward  ends  of  the  chamber,  represent  small 
differences  between  large  numbers  and  are  not  appropriately  compared  in  the 
same  fashion.  Rather,  we  choose  to  compare  the  entire  pressure-difference 
versus  time  profiles  in  Figure  6.  Similarly,  we  display  a  comparison  of  axial 
flamespread  profiles  in  Figure  7.  We  note  a  favorable  level  of  agreement  for 
all  parameters. 


TABLE  1.  SUMMARY  OF  NOVA /TDNOVA  RESULTS 


CODE 


MAX  PRESSURE 
(MPa) 


MUZZLE  VELOCITY 
(m/s) 


INITIAL  REVERSE 
PRESSURE  DIFF 
(MPa) 


NOVA  (30x1  mesh) 

323 

814 

-26 

NOVA  -  modified 

treatment  of  end 
ullage  (30x1  mesh) 

325 

816 

-21 

TDNOVA  (30x7  mesh) 

326 

820 

-21 

TDNOVA  (30x1  mesh) 

327 

820 

-21 
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Figure  4.  Quasi-Two-Dimensional  Mesh  of  TDNOVA 
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Figure  5.  Schematic  Representation  of  Quasi-One-Dimensional  Test  Problem 
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Figure  6.  NOVA  and  TDNOVA  Predictions  of  Pressure  Difference  Versus  Time 
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c.  Influence  of  Criterion 

We  noted  earlier  that,  upon  completion  of  flamespread,  rupture  of  any  bag 
sidewall  material  present,  and  equilibration  of  the  radial  structure  of  the 
pressure  field  to  within  some  user-specified  limit,  TDNOVA  introduces  a  quasi- 
two-dimensional  approach  for  the  duration  of  the  calculation.  This 
transformation  is  invoked  in  the  interest  of  economy  and,  in  consideration  of 
the  scope  for  which  TDNOVA  is  intended,  appears  to  be  well  motivated.  The 
particular  criterion  employed  involves  Ptol,  a  parameter  which  is  compared,  at 
each  axial  location,  to  the  difference  in  values  for  the  pressure  at  the  tube 
wall  and  centerline,  divided  by  the  value  at  the  tube  wall.  When  this  # 

quantity  becomes  less  than  at  all  axial  stations,  the  transformation 

takes  place. 

As  this  transformation  carries  with  it  a  number  of  assumptions^  required 
to  establish  a  quasi-one-dimensional  description  of  flow  within  the  two-phase 
medium,  we  undertook  to  determine  the  influence  of  Pto^  on  the  remainder  of 
the  solution.  A  previously  established  data  base4  (see  Figure  8  and  Appendix 
B)  for  the  155-mm,  M203  Propelling  Charge  was  employed,  with  Pto^  varied 
over  a  wide  range  of  values.  Results  are  summarized  in  Table  2,  and  a 
comparison  of  pressure-difference  versus  time  profiles  is  displayed  in  Figure 
9.  All  values  of  Ptol  GSual  to  or  greater  than  0.05  yielded  identical 
results,  as  the  time  for  the  final  point  of  bag  rupture  became  the  controlling 
parameter,  prohibiting  transformation  at  any  earlier  times.  Even  for  values 
Ptol  as  smaH  as  0.005,  results  remain  virtually  unchanged.  Pressure  and 
temperature  fields,  at  the  instant  of  transformation,  are  compared  in  Figures 
10  and  11.  Little  difference  in  structure  is  noted  for  the  various 
conditions.  Prediction  of  flamespread  is,  of  course,  totally  unaffected  by 
PtoP  completion  being  a  requirement  for  transformation. 


TABLE  2.  INFLUENCE  OF  Ptol  ON  TDNOVA  RESULTS 


Ptol  (and 
time  of  trans¬ 
formation,  ms) 

MAX  PRESSURE 
(MPa) 

MUZZLE  VELOCITY 
(m/s) 

INITIAL  REVERSE 
PRESSURE  DIFF 

(MPa) 

0.005 

(3.51) 

365 

842 

-3 

0.010 

(3.49) 

363 

842 

-3 

0.050 

(3.19) 

363 

842 

-2 

>0.050 

(3.19) 

363 

842 

-2 
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Figure  8.  Schematic  Representation  of  155-mm,  M203  Propelling  Charge 


Figure  9.  Influence  of  P-^ol  on  TDNOVA  Predictions  of  Pressure  Difference 
Versus  Time 
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Figure  10.  TDNOVA  Predictions  of  Pressure  Fields  at  Times  of  Transformation 
for  Several  Values  of  P^ol 
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Figure  11.  TDNOVA  Predictions  of  Temperature  Fields  at  Times  of 
Transformation  for  Several  Values  of  Ptol 
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D.  Influence  of.  Mesh 


To  complete  this  baseline  evaluation  of  TDNOVA,  attention  must  also  be 
directed  to  the  sensitivity  of  solutions  to  the  mesh  employed  to  represent  the 
two-dimensional,  two-phase  region  of  flow.  The  input  data  base  provided  in 
Appendix  B  was  again  selected  for  this  phase  of  the  study.  From  10  to  35 
axial  mesh  points  and  from  3  to  9  radial  mesh  points  were  employed  to 
represent  the  region  occupied  by  the  propellant  bed.  While  not  fundamentally 
linked  to  the  limitations  imposed  by  the  macroscopic  nature  of  the  governing 
equations  for  TDNOVA,  a  selection  of  mesh  size  somewhere  in  the  range  studied 
is  certainly  compatible  with  both  the  intended  purpose  and  physical  scope  of 
TDNOVA.  Further,  in  the  case  of  any  extensive  propelling  charge  design 
studies,  the  need  for  economy  may  also  limit  one  to  this  range  of  values. 

A  summary  of  results  from  these  calculations  is  provided  in  Table  3.  We 
note  some  apparent  dependence  of  predicted  performance  on  the  number  of  radial 
mesh  points,  though  the  total  spread  is  less  than  45s  for  values  of  maximum 
chamber  pressure  and  less  than  2%  for  those  of  muzzle  velocity.  As  before, 
the  initial  reverse  pressure  difference,  being  calculated  as  the  difference  of 
two  large  numbers,  exhibits  a  large  percentage  but  small  absolute  variation. 
Selected  pressure-difference  versus  time  profiles,  flamespread  contours,  and 
pressure  field  plots  are  displayed  in  Figures  12  through  16. 


TABLE  3.  SUMMARY  OF  TDNOVA  RESULTS  FOR  VARIOUS  MESHES 


MESH 

MAX  PRESSURE 

MUZZLE  VELOCITY 

INITIAL  REVERSE 

axial  pts  x 

(MPa) 

(m/s) 

PRESSURE  DIFF 

radial  pts) 

(MPa) 

20x3 

357 

836 

-3 

25x3 

354 

834 

~3 

30x3 

357 

836 

-3 

35x3 

358 

836 

-4 

10x5 

359 

838 

-3 

20x5 

359 

838 

-4 

25x5 

363 

839 

~3 

30x5 

359 

830 

-5 

35x5 

356 

831 

-7 

20x7 

362 

840 

-2 

25x7 

364 

843 

-2 

30x7 

363 

842 

“3 

35x7 

366 

835 

-3 

20x9 

364 

842 

-2 

25x9 

368 

845 

-1 
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Figure  12.  TDNOVA  Predictions  of  Pressure  Difference  Versus  Time  for 
Various  Meshes 


RADIAL  LOCATION  (cm) 
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20x5  MESH 


Note:  flame  contours  displayed  every  0.24  ms 


Figure  13.  TDNOVA  Predictions  of  Flamespread  for  Various  Meshes 
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Figure  14.  TDNOVA  Prediction  of  Pressure  Fields  for  a  20x5  Mesh 
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Figure  15.  TDNOVA  Prediction  of  Pressure  Fields  for  a  25x9  Mesh 
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Figure  16.  TDNOVA  Prediction  of  Pressure  Fields  for  a  35x7  Mesh 


III.  CONCLUSIONS 


In  an  effort  to  provide  a  baseline  evaluation  of  the  TDNOVA  code,  several 
series  of  calculations  were  performed.  Based  on  the  results  of  these 
calculations,  the  following  conclusions  can  be  drawn: 

1.  Simulations  of  a  quasi-one-dimensional  propelling  charge  obtained 
using  TDNOVA  and  NOVA,  its  quasi-one-dimensional  predecessor,  are  essentially 
equivalent.  When  TDNOVA  is  modified  to  introduce  an  immediate  transformation 
to  its  quasi-two-dimensional  mode  and  when  NOVA  is  also  modified  to  maintain  a 
lumped-parameter  representation  of  regions  of  axial  ullage,  results  from  the 
two  codes  become  virtually  identical. 

2.  For  at  least  one,  relevant,  bagged-charge  problem,  results  provided 
by  TDNOVA  are  only  minimally  influenced  by  the  value  selected  for  Ptoi»  a 
parameter  used  to  identify  adequate  equilibration  of  the  radial  pressure  field 
before  transformation  to  a  quasi-two-dimensional  representation  of  flow  is 
allowed.  Antecedent  requirements  for  bag  rupture  and  completion  of 
f lamespread  apparently  allow  substantial  equilibration  of  radial  pressures 
prior  to  application  of  the  criterion. 

3.  A  limited  study  of  the  influence  of  mesh  density  on  TDNOVA  results 
failed  to  demonstrate  absolute  convergence  of  results;  nevertheless,  the 
variation  in  predicted  quantities  was  shown  to  be  acceptably  small  for  a 
number  of  meshes  covering  the  current  range  of  practical  interest. 
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APPENDIX  A 

INPUT  DATA  FOR  NOVA /TDNOVA  COMPARISON  CALCULATIONS 
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NOTE:  The  NOVA  code  was  developed  prior  to  the  introduction  of 
the  International  System  of  Units  (SI)  and  employs  English  units 
throughout.  A  summary  of  conversion  factors  required  to  effect 
conversion  to  SI  units,  employed  in  TDNOVA ,  is  provided  below. 


TO  CONVERT 

FROM  NOVA 

UNITS 

TO  TDNOVA 
UNITS 

MULTIPLY  BY 

LENGTH 

in  . 

cm 

2 . 54 

MASS 

lbm 

gm 

453 . 59237 

TEMPERATURE 

°R 

K 

5/9 

FORCE 

lbf 

N 

4 . 448222 

VELOCITY 

in . /sec 

cm/ sec 

2 . 54 

PRESSURE 

lbf /in . 2 

MPa 

0 . 006894757 

DENSITY 

lbm/ in . 2 

gm/ cm3 

27 . 679905 

COVOLUME 

in . 2 /lbm 

cm3//9m 

0 .036127292 

INTERNAL  ENERGY 

lbf -in . / lbm 

J/gm 

0 . 000249089 

THERMAL  CONDUCTIVITY 

lbf -in . /in . 

-sec-°R  J/cm- 

sec-K  0.080068 

THERMAL  DIFFUSIVITY 

in . 2/sec 

cm  /sec 

6 . 4516 

BURN  RATE  PRE -EXPONENT  in./sec-psi 

n  cm/sec-MPa 

n  _ 2^.54 _ 

(0.00689457) 
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MQJZA  CQDEi  CQMEABISQM  IMEUI  DATA  EASE 


CONTROL  PARAMETERS 

PRINT  T 

GRAPH  T 

DISK  WRITE  F 

DISK  READ  F 

I.B.  TABLE  T 

FLAME  TABLE  T 

PRESSURE  TABLES  T 

EROSIVE  EFFECT  0 

DYNAMIC  EFFECT  0 

WALL  TEMPERATURE  CALCULATION  0 

LEFT  HAND  BOUNDARY  CONDITION  0 

RIGHT  HAND  BOUNDARY  CONDITION  0 

LEFT  HAND  RESERVOIR  0 

RIGHT  HAND  RESERVOIR  0 

BED  PRECOMPRESSED  0 

HEAT  LOSS  CALCULATION  0 

INSULATING  LAYER  0 

BORE  RESISTANCE  FUNCTION  0 

EXPLICIT  COMPACTION  WAVE  0 

MUZZLE  BLOWDOWN  ANALYSIS  0 

CALCOMP  SUMMARY  PLOTS  0 

INTEGRATION  PARAMETERS 

NUMBER  OF  STATIONS  AT  WHICH  DATA  ARE  STORED  30 

NUMBER  OF  STEPS  BEFORE  LOGOUT  100 

TIME  STEP  FOR  DISK  START  0 

NUMBER  OF  STEPS  FOR  TERMINATION  3000 

TIME  BEFORE  PRINTOUT  .0005 

PRESSURE  RATIO  FOR  LP  ANALYSIS  OF  LARGE 

ULLAGE  REGION  . 2 

TIME  FOR  TERMINATION  (SEC)  .05 

PROJECTILE  TRAVEL  FOR  TERMINATION  (INS)  205. 

MAXIMUM  TIME  STEP  (SEC)  .0001 

STABILITY  SAFETY  FACTOR  2. 

SOURCE  STABILITY  FACTOR  .05 

SPATIAL  RESOLUTION  FACTOR  .01 

TIME  INTERVAL  FOR  I.B.  TABLE  STORAGE  (SEC)  .0002 

TIME  INTERVAL  FOR  PRESSURE  TABLE  STORAGE  (SEC)  .0002 

FILE  COUNTERS 

NUMBER  OF  STATIONS  TO  SPECIFY  TUBE  RADIUS  3 

NUMBER  OF  TIMES  TO  SPECIFY  PRIMER  DISCHARGE  3 

NUMBER  OF  POSITIONS  TO  SPECIFY  PRIMER  DISCHARGE  3 
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FILE  COUNTERS  (CONTINUED) 
NUMBER  OF  ENTRIES  IN  BORE  RESISTANCE  TABLE 
NUMBER  OF  ENTRIES  IN  WALL  TEMPERATURE  TABLE 
NUMBER  OF  ENTRIES  IN  FILLER  ELEMENT  TABLE 
NUMBER  OF  TYPES  OF  PROPELLANTS 
NUMBER  OF  BURN  RATE  DATA  SETS 
NUMBER  OF  ENTRIES  IN  VOID  FRACTION  TABLES 
NUMBER  OF  ENTRIES  IN  PRESSURE  HISTORY  TABLES 
NUMBER  OF  ENTRIES  IN  LEFT  BOUNDARY  SOURCE  TABLE 
NUMBER  OF  ENTRIES  IN  RIGHT  BOUNDARY  SOURCE  TABLE 
NUMBER  OF  WALL  STATIONS  FOR  INVARIANT  EMBEDDING 
NUMBER  OF  BED  STATIONS  FOR  INVARIANT  EMBEDDING 


7 

0 

0 

1 

2 

0 

3 

0 

0 

0 

0 


0 


0 


GENERAL  PROPERTIES  OF  INITIAL  AMBIENT  GAS 


INITIAL  TEMPERATURE  (DEG  R)  530. 

INITIAL  PRESSURE  (PSI)  14.7 

MOLECULAR  WEIGHT  ( LBM/LBMOL )  23.36 

RATIO  OF  SPECIFIC  HEATS  1.243 


GENERAL  PROPERTIES  OF  PROPELLANT  BED 

INITIAL  TEMPERATURE  (DEG  R)  530. 

VIRTUAL  MASS  COEFFICIENT  FOR  MOMENTUM  TRANSFER  0. 

VIRTUAL  MASS  COEFFICIENT  FOR  ENERGY  DISSIPATION  0. 

MINIMUM  IMPACT  VELOCITY  FOR  EXPLICIT 

COMPACTION  WAVE  (IN/SEC)  100000000. 

FRICTION  FACTOR  1- 


PROPERTIES  OF  PROPELLANT  1 


PROPELLANT  TYPE 
MASS  OF  PROPELLANT  (LBM) 
DENSITY  OF  PROPELLANT  (LBM) 
FORM  FUNCTION  INDICATOR 
OUTSIDE  DIAMETER  (INS) 
INSIDE  DIAMETER  (INS) 

LENGTH  (INS) 

NUMBER  OF  PERFORATIONS 
SLOT  WIDTH  (INS) 


M30A1 ,  RAD - 7 7G -06 9 80 5 
26  .  15 
.  0572 
7 

.4173 
.0338 
.  9481 
7  . 

0  . 


RHEOLOGICAL  PROPERTIES 

SPEED  OF  COMPRESSION  WAVE  IN  SETTLED  BED 

(IN/SEC)  6000. 

SETTLING  POROSITY 

SPEED  OF  EXPANSION  WAVE  (IN/SEC)  50000. 
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SOLID  PHASE  THERMOCHEMISTRY 


MAXIMUM  PRESSURE  FOR  BURN  RATE  DATA 
( LBF/IN**2 ) 

BURNING  RATE  PRE-EXPONENTIAL  FACTOR 
( IN/SEC/PSI**BN) 

BURNING  RATE  EXPONENT 
MAXIMUM  PRESSURE  FOR  BURN  RATE  DATA 
( LBF / IN*  *  2 ) 

BURNING  RATE  PRE -EXPONENTIAL  FACTOR 
( IN/SEC/PSI **BN ) 

BURNING  RATE  EXPONENT 
BURNING  RATE  CONSTANT 
IGNITION  TEMPERATURE  (DEG  R) 

ARRHENIUS  ACTIVATION  ENERGY  ( LBF -IN/LBMOL ) 
FREQUENCY  FACTOR  (SEC**-1) 

THERMAL  CONDUCTIVITY  ( LBF /SEC /DEG  R) 
THERMAL  DIFFUSIVITY  <IN**2/SEC) 

EMMIS IVITY  FACTOR 


lOOOO . 


.  006918 
.  6337 


100000 . 

.001743 
.  7864 
0  . 

800  . 

0  . 

0  . 

.02 
.  0003 
0  . 


GAS  PHASE  THERMOCHEMISTRY 

CHEMICAL  ENERGY  RELEASED  IN  BURNING 

( LBF-IN/LBM )  17600000. 

MOLECULAR  WEIGHT  ( LBM/LBMOL )  23.36 

RATIO  OF  SPECIFIC  HEATS  1 • 24 

COVOLUME  28.50 


PACKAGE 

1 


LOCATION  OF  PACKAGE ( S ) 
LEFT  BDDY  (INS)  RIGHT  BDDY  (INS) 
1.00  31.00 


MASS  ( LBM ) 
26.15 


PROPERTIES  OF  IGNITER 

CHEMICAL  ENERGY  RELEASED  IN  BURNING 

(LBF-IN/LBM)  10000000. 

MOLECULAR  WEIGHT  (LBM/LBMOL)  23.36 

RATIO  OF  SPECIFIC  HEATS  1 • 24 

SPECIFIC  VOLUME  OF  SOLID  (IN**3/LBM)  0. 


IGNITER  DISCHARGE  FUNCTION  (LBM/IN/SEC) 


POS.  (INS) 

0 . 00 

0 . 98 

0 . 99 

TIME  (SEC) 
0.000 

6 . 00 

6 . 00 

0 . 00 

0 .010 

6 . 00 

6 . 00 

0.00 

0.011 

0 . 00 

0 . 00 

0 . 00 
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PARAMETERS  TO  SPECIFY  TUBE  GEOMETRY 


DISTANCE  (INS) 

RADIUS  (INS) 

0 . 00 

3 . 50 

35 . 00 

3 .09 

240 . 00 

3 .09 

BORE  RESISTANCE 

POSITION  (INS) 

RESISTANCE  (PSD 

35 . 00 

250  . 

35 . 40 

3350  . 

36 . 00 

4950  . 

36 . 55 

3625  . 

37.05 

3250  . 

39 . 50 

2500  . 

240 . 00 

1500  . 

TABLE 


THERMAL  PROPERTIES  OF  TUBE 


THERMAL  CONDUCTIVITY  (LBF/SEC/DEG  R)  0. 

THERMAL  DIFFUSIVITY  (IN**2/SEC)  0. 

EMISSIVITY  FACTOR  0. 

INITIAL  TEMPERATURE  (DEG  R)  530. 

PROJECTILE  AND  RIFLING  DATA 

INITIAL  POSITION  OF  BASE  OF  PROJECTILE  (INS)  35.00 

MASS  OF  PROJECTILE  ( LBM )  103.00 

POLAR  MOMENT  OF  INERTIA  (LBM-IN**2)  0. 

ANGLE  OF  RIFLING  (DEG)  0. 


POSITIONS  FOR  PRESSURE  TABLE  STORAGE  (INS) 
0.10  17.50  34.75 
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IDtLQJZA  QQUEj.  CQMEABISQN  IMEUI  DAIA  BASE 


CONTROL  PARAMETERS 


NPRINT  ( 0  =NO  PRINT, 1=PRINT)  1 
NSUMRY  ( 0  =NO  SUMMARY  TABLES , 1 = YES )  1 
NPLOT  (0=NO  ISOMETRIC  CARPET  PLOTS , 1 =PLOT )  1 
NVHL  ( 0=HIDDEN  LINES  DELETED , 1 =RETAINED )  0 
NPLCON  ( 0=N0  CONTOUR  PLOTS , 1 =PLOT )  0 
NPLFLO  ( 0=NO  FLOW  PLOTS , 1 =PLOT )  0 
NPLFLM  ( 0=NO  FLAMESPREAD  PLOT,l=PLOT)  1 
NDSKW  ( 0=NO  DISC  SAVE,1=DISC  SAVE)  0 
NDSKR  ( 0=N0  DISC  START , >0=DISC  START  AT 

STEP  NDSKR)  0 


ISOMETRICALLY  PLOTTED  QUANTITIES  <0=NO,1=YES) 


MESH  0 
POROSITY  0 
GRANULAR  STRESS  0 
PRESSURE  1 
DENSITY  0 
GAS  AXIAL  VELOCITY  0 
SOLID  AXIAL  VELOCITY  0 
GAS  RADIAL  VELOCITY  0 
SOLID  RADIAL  VELOCITY  0 
GAS  TEMPERATURE  1 
PARTICLE  SURFACE  TEMPERATURE  0 


CONTOUR  PLOTTED  QUANTITIES  <0=N0,1=YES) 


MESH  0 

POROSITY  0 

GRANULAR  STRESS  0 

PRESSURE  0 

DENSITY  0 

GAS  AXIAL  VELOCITY  0 

SOLID  AXIAL  VELOCITY  0 

GAS  RADIAL  VELOCITY  0 

SOLID  RADIAL  VELOCITY  0 

GAS  TEMPERATURE  0 

PARTICLE  SURFACE  TEMPERATURE  0 

SCALE  FACTOR  FOR  PLOTTING  .4 

LENGTH  OF  Z-AXIS  IN  CALCOMP  PLOTS  (INS)  12. 

LENGTH  OF  R-AXIS  (INS)  4. 

LENGTH  OF  ORDINATE  AXIS  (INS)  5. 
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LOGOUT  PARAMETERS 


NUMBER  OF  STEPS  BEFORE  LOGOUT  2000 

TIME  INCREMENT  BEFORE  LOGOUT  (MSEC)  .5 

NUMBER  OF  PRESSURE  SUMMARY  STATIONS  3 

TIME  INCREMENT  FOR  PRESSURE  SUMMARY  STORAGE 

(MSEC)  .2 

TERMINATION  PARAMETERS 

MAXIMUM  NUMBER  OF  STEPS  BEFORE  TERMINATION  2000 

MAXIMUM  INTEGRATION  TIME  (MSEC)  25. 

MAXIMUM  PROJECTILE  TRAVEL  (CMS)  520.7 

MESH  PARAMETERS 

MESH  ALLOCATION  MODE  ( 0=STATIC , 1 = DYNAMIC )  0 

MAXIMUM  NUMBER  OF  STORAGE  POINTS  FOR  DYNAMIC 

MESH  ALLOCATION  0 

NUMBER  OF  MESH  POINTS  IN  AXIAL  DIRECTION  30 

NUMBER  OF  MESH  POINTS  IN  RADIAL  DIRECTION  7 

NUMBER  OF  ITERATIONS  TO  DETERMINE  INITIAL 

MESH  200 

SAFETY  FACTOR  FOR  C-F-L  CRITERION  1.1 

MAXIMUM  FRACTIONAL  DISPLACEMENT  FOR  CONVER¬ 
GENCE  OF  INITIAL  MESH  DISTRIBUTION  .00001 

OVER-RELAXATION  FACTOR  FOR  DETERMINATION  OF 

INITIAL  MESH  DISTRIBUTION  1.6 

PRESSURE  TOLERANCE  FACTOR  FOR  REDUCTION  TO 

QUASI -TWO-DIMENSIONAL  REPRESENTATION  .05 

AXIAL  SPATIAL  RESOLUTION  FACTOR  .1 

RADIAL  SPATIAL  RESOLUTION  FACTOR 

AMBIENT  CONDITIONS 

INITIAL  TEMPERATURE  (DEG  K)  294.4 

INITIAL  PRESSURE  (MPA)  !l014 

CHARGE  STANDOFF  (CMS)  0! 

SOLID  PHASE  CONSTITUTIVE  DATA 

INITIAL  MASS  OF  GRANULAR  BED  (KG)  11  86 

INITIAL  POROSITY  OF  GRANULAR  BED  0 

SETTLING  POROSITY  OF  GRANULAR  BED  4243 

SPEED  OF  COMPRESSION  WAVE  (M/SEC)  152  4 

SPEED  OF  EXPANSION  WAVE  (M/SEC)  1270  0 

DENSITY  OF  SOLID  PHASE  ( GM/CC )  1  583 

THERMAL  CONDUCTIVITY  ( J/CM-SEC -DEG  K)  0016 

THERMAL  DIFFUSIVITY  (CM**2/SEC)  0006 
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GAS  PHASE  CONSTITUTIVE  DATA 


RATIO  OF  SPECIFIC  HEATS  1.243 

MOLECULAR  WEIGHT  (GM/GM-MOL)  23.36 

COVOLUME  (CC/GM)  1.030 


SOLID  PHASE  COMBUSTION  CHARACTERISTICS 

IGNITION  TEMPERATURE  (DEG  K)  444.4 

CHEMICAL  ENERGY  (J/GM)  4384. 


MAX  PRESSURE  ADD  CONSTANT  PRE-EXPONENT  EXPONENT 


(MPA) 

(CM/SEC) 

<CM/SEC-MPA**BN) 

68 . 95 

0  . 

.4117 

.  6337 

689 . 50 

0  . 

.  2218 

.  7864 

EXTERNAL  DIAMETER  (CM) 
LENGTH  (CM) 

GRAIN  GEOMETRY 

1 . 060 
2 . 408 

DIAMETER  OF  PERFORATIONS 
NUMBER  OF  PERFORATIONS 

(CM) 

.  086 

7  . 

AX  POS 

RAD  POS 

FLOW  RES 

REACTIVITY 

NO.  PTS 

DATA  TYPE 

(CMS) 

(CMS  ) 

DATA 

DATA 

PRE-ASS . 

( 0  =  D , 1 =N ) 

CONFIGURATION 

OF 

REAR  OF 

BAG 

2. 54 

0 . 00 

0 

0 

0 

0 

2. 54 

8.86 

0 

0 

0 

0 

CONFIGURATION 

OF 

FRONT  OF 

BAG 

78. 74 

0 . 00 

0 

0 

0 

0 

78 . 74 

7 . 97 

0 

0 

0 

0 

CONFIGURATION 

OF 

INSIDE  OF 

BAG 

2 . 54 

0.00 

0 

0 

0 

0 

78. 74 

0.00 

0 

0 

0 

0 

CONFIGURATION 

OF 

OUTSIDE  OF  BAG 

2 . 54 

8. 86 

0 

0 

0 

0 

78 . 74 

7 . 97 

0 

0 

0 

0 
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AXIAL  POSITION 
<  CMS  ) 


RADIAL  POSITION 
(CMS  ) 


CONFIGURATION  OF  BREECH 

0.00  0.00 

0.00  8.89 


CONFIGURATION  OF  PROJECTILE  BASE 

0 . 00 
7 . 85 


CONFIGURATION  OF  INSIDE  BOUNDARY 

0 . 00 

0.00 


CONFIGURATION  OF  OUTSIDE  BOUNDARY 

8 . 89 
7 . 85 


REPRESENTATION  OF  IGNITION  TRAIN 

NCCORE  ( 0  =  NO  CENTERCORE , 1=YES ) 

BASEPAD  REACTIVITY  DATA 
NTABAG  ( 0  =NO  EXTERNAL  STIMULUS , 1 =YES )  1 


CHEMICAL  ENERGY  OF  EXTERNAL  STIMULUS  (J/GM)  2491. 


0 . 00 
88 . 90 


0 . 00 
88 . 90 


88 . 90 
88 . 90 


RATE  OF  DISCHARGE  OF  EXTERNAL  STIMULUS 


RATE  (GM/CC/SEC)  AT  TIME  0.00  MSEC 


AXIAL  LOCATION 
0 . 00 
2 . 50 
2 . 52 


RADIAL 
0 . 00 

(  CM) 

4  .  65 
4  .  65 
0 . 00 


LOCATION  (CM) 


8  . 

82 

8  . 

.  84 

4  . 

.  65 

0  . 

.  00 

4  . 

.  65 

0  . 

.  00 

0  . 

.  00 

0  . 

.  00 
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o  o 


RATE  OF  DISCHARGE  OF  EXTERNAL  STIMULUS  (CONTINUED) 
RATE  (GM/CC/SEC)  AT  TIME  10.00  MSEC 

RADIAL  LOCATION  (CM) 


0.00 

8 . 82 

00 

00 

AXIAL 

LOCATION  (CM) 

0.00 

4 . 65 

4 . 65 

0.00 

2 . 50 

4 . 65 

4 . 65 

0.00 

2. 52 

0.00 

0 . 00 

0.00 

RATE  ( 

GM/CC/SEC)  AT 

TIME  11. 

00  MSEC 

RADIAL 

LOCATION 

(CM) 

0 . 00 

8.82 

8 . 84 

AXIAL 

LOCATION  (CM) 

0.00 

0.00 

0.00 

0 . 00 

2.50 

0.00 

0.00 

0.00 

2.52 

0 . 00 

0.00 

0.00 

PROPERTIES  OF  PROJECTILE 

PROJECTILE  MASS  (KG) 

number  of  entries  in  bore  resistance  table 

RESISTANCE  LAM  NUMBER 
NUMBER  OF  FILLER  ELEMENTS 


BORE  RESISTANCE  DATA 


PROJECTILE  TRAVEL 
(CMS  ) 

0.000 
1.016 
2.540 
3. 937 
5 . 207 
11.430 
520 . 700 


RESISTIVE  PRESSURE 
(MPA) 

1 . 72 

23.10 

34 . 10 
25 . 00 
22 . 40 
17.20 
10. 30 


LOCATION  OF  POINTS  FOR  PRESSURE 

AXIAL  LOCATION  HALL  (0)  OR  AXIS 

(  CMS  ) 


SUMMARY  TABLE 
(  1  ) 


.  25 
44 . 45 
88 . 27 
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0 

0 

0 


APPENDIX  B 

INPUT  DATA  FOR  TDNOVA  Ptol  AND  MESH-SENSITIVITY  CALCULATIONS 
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IDNQJZA  CQDEjl  Etfii  AMD  MESH  SIUDI  DAIA  BASE 


CONTROL  PARAMETERS 


NPRINT  (0=NO  PRINT, 1=PRINT)  1 
NSUMRY  ( 0  =NO  SUMMARY  TABLES , 1 =YES )  2 
NPLOT  ( 0=N0  ISOMETRIC  CARPET  PLOTS , 1 =PLOT )  1 
NVHL  ( 0=HIDDEN  LINES  DELETED , 1 =RETAINED )  0 
NPLCON  ( 0  =  NO  CONTOUR  PLOTS , 1 =PLOT )  0 
NPLFLO  ( 0  =NO  FLOW  PLOTS , 1 =PLOT )  0 
NPLFLM  ( 0  =  NO  FLAMESPREAD  PLOT,l=PLOT)  1 
NDSKW  ( 0=NO  DISC  SAVE,1=DISC  SAVE)  0 
NDSKR  ( 0  =  NO  DISC  START , >0  =  DISC  START  AT 

STEP  NDSKR)  0 


ISOMETRICALLY  PLOTTED  QUANTITIES  (0=NO,1=YES) 


MESH  0 

POROSITY  0 

GRANULAR  STRESS  0 

PRESSURE  1 

DENSITY  0 

GAS  AXIAL  VELOCITY  0 

SOLID  AXIAL  VELOCITY  0 

GAS  RADIAL  VELOCITY  0 

SOLID  RADIAL  VELOCITY  0 

GAS  TEMPERATURE  1 

PARTICLE  SURFACE  TEMPERATURE  0 

CONTOUR  PLOTTED  QUANTITIES  (0=NO,1=YES) 

MESH  0 

POROSITY  0 

GRANULAR  STRESS  0 

PRESSURE  0 

DENSITY  0 

GAS  AXIAL  VELOCITY  0 

SOLID  AXIAL  VELOCITY  0 

GAS  RADIAL  VELOCITY  0 

SOLID  RADIAL  VELOCITY  0 

GAS  TEMPERATURE  0 

PARTICLE  SURFACE  TEMPERATURE  0 

SCALE  FACTOR  FOR  PLOTTING  .4 

LENGTH  OF  Z-AXIS  IN  CALCOMP  PLOTS  (INS)  12. 

LENGTH  OF  R-AXIS  (INS)  4. 

LENGTH  OF  ORDINATE  AXIS  (INS)  5. 
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LOGOUT  PARAMETERS 


NUMBER  OF  STEPS  BEFORE  LOGOUT  2000 

TIME  INCREMENT  BEFORE  LOGOUT  (MSEC)  .5 

NUMBER  OF  PRESSURE  SUMMARY  STATIONS  3 

TIME  INCREMENT  FOR  PRESSURE  SUMMARY  STORAGE 

(MSEC)  .2 


TERMINATION  PARAMETERS 

MAXIMUM  NUMBER  OF  STEPS  BEFORE  TERMINATION  2000 

MAXIMUM  INTEGRATION  TIME  (MSEC)  25. 

MAXIMUM  PROJECTILE  TRAVEL  (CMS)  520.7 


MESH  PARAMETERS 


MESH  ALLOCATION  MODE  ( 0=STATIC , 1 =DYNAMIC )  0 

MAXIMUM  NUMBER  OF  STORAGE  POINTS  FOR  DYNAMIC 

MESH  ALLOCATION  0 

NUMBER  OF  MESH  POINTS  IN  AXIAL  DIRECTION  VARIOUS 

NUMBER  OF  MESH  POINTS  IN  RADIAL  DIRECTION  VARIOUS 

NUMBER  OF  ITERATIONS  TO  DETERMINE  INITIAL 

MESH  200 

SAFETY  FACTOR  FOR  C-F-L  CRITERION  1.1 

MAXIMUM  FRACTIONAL  DISPLACEMENT  FOR  CONVER¬ 
GENCE  OF  INITIAL  MESH  DISTRIBUTION  .00001 

OVER-RELAXATION  FACTOR  FOR  DETERMINATION  OF 

INITIAL  MESH  DISTRIBUTION  1.6 

PRESSURE  TOLERANCE  FACTOR  FOR  REDUCTION  TO 

QUASI-TWO-DIMENSIONAL  REPRESENTATION  VARIOUS 

AXIAL  SPATIAL  RESOLUTION  FACTOR  .1 

RADIAL  SPATIAL  RESOLUTION  FACTOR  .1 

AMBIENT  CONDITIONS 

INITIAL  TEMPERATURE  (DEG  K)  294.4 

INITIAL  PRESSURE  (MPA)  .1014 

CHARGE  STANDOFF  (CMS)  0. 


SOLID  PHASE  CONSTITUTIVE  DATA 


INITIAL  MASS  OF  GRANULAR  BED  (KG)  11.86 

INITIAL  POROSITY  OF  GRANULAR  BED  0. 

SETTLING  POROSITY  OF  GRANULAR  BED  0. 

SPEED  OF  COMPRESSION  WAVE  (M/SEC)  152.4 

SPEED  OF  EXPANSION  WAVE  (M/SEC)  1270.0 

DENSITY  OF  SOLID  PHASE  (GM/CC)  1.583 

THERMAL  CONDUCTIVITY  ( J/CM-SEC -DEG  K)  .0016 

THERMAL  DIFFUSIVITY  <CM**2/SEC)  .0006 
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GAS  PHASE  CONSTITUTIVE  DATA 


RATIO  OF  SPECIFIC  HEATS 
MOLECULAR  WEIGHT  (GM/GM-MOL) 
COVOLUME  (CC/GM) 


SOLID  PHASE  COMBUSTION  CHARACTERISTICS 


IGNITION  TEMPERATURE  (DEG  K) 
CHEMICAL  ENERGY  (J/GM) 


MAX  PRESSURE 
(MPA) 


ADD  CONSTANT 
(CM/SEC) 


PRE-EXPONENT  EXPONENT 

(CM/SEC-MPA**BN) 


68.95  0. 

689.50  0. 


•4117  .6337 

.2218  .7864 


GRAIN  GEOMETRY 


EXTERNAL  DIAMETER  (CM) 

LENGTH  (CM) 

DIAMETER  OF  PERFORATIONS  (CM) 
NUMBER  OF  PERFORATIONS 


1.060 
2 . 408 
.086 
7. 


AX  POS 
(CMS  ) 


2.54 

2.54 

2.54 


RAD  POS 

FLOW  RES 

REACTIVITY 

NO.  PTS 

DATA  TYPE 

(CMS  ) 

DATA 

DATA 

PRE-ASS . 

(0=D, 1=N) 

CONFIGURATION 

OF  REAR  OF 

BAG 

1 . 27 

1 

1 

0 

0 

3.40 

1 

0 

0 

0 

7 . 62 

0 

0 

0 

0 

CONFIGURATION  OF  FRONT  OF  BAG 

78.74  1.27  2  0 

78.74  7.62  0  0 


0 

0 


CONFIGURATION  OF  INSIDE  OF  BAG 

2.54  1.27  3  0 

78.74  1.27  0  0 


0 

0 


2.54 
27 . 94 
78 . 74 


CONFIGURATION 

7 . 62 

7.62 

7.62 


OF  OUTSIDE  OF  BAG 

0 
0 
0 


0 

0 

0 


1 

4 

0 
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0 

0 

0 


AXIAL  POSITION 
(CMS) 


RADIAL  POSITION 
(CMS  ) 


CONFIGURATION  OF  BREECH 

0 . 00 
0 . 00 
-3.  45 


0 . 00 
7 .06 
8 . 48 


CONFIGURATION  OF  PROJECTILE  BASE 

87  .  38 
87. 38 
96 . 42 


0 . 00 
7.14 
7 . 85 


CONFIGURATION  OF  INSIDE  BOUNDARY 

0 . 00 
87  .  38 


0.00 
0 . 00 


CONFIGURATION  OF  OUTSIDE  BOUNDARY 

-3. 45 
92 . 46 
96. 42 


8 . 48 
8.05 
7 . 85 


REPRESENTATION  OF  IGNITION  TRAIN 

NCCORE  ( 0  =  N0  CENTERCORE , 1=YES )  1 
BASEPAD  REACTIVITY  DATA  1 
NTABAG  ( 0=NO  EXTERNAL  STI MULUS , 1 = YES )  0 


PROPERTIES  OF  SOLID  PHASE  IN  CENTERCORE 


SOLID  PHASE  CONSTITUTIVE  DATA 


INITIAL  MASS  OF  GRANULAR  BED  (KG)  .1134 

INITIAL  POROSITY  OF  GRANULAR  BED  0. 

SETTLING  POROSITY  OF  GRANULAR  BED  .40 

SPEED  OF  COMPRESSION  WAVE  (M/SEC)  442. 

SPEED  OF  EXPANSION  WAVE  (M/SEC)  1270. 

DENSITY  OF  SOLID  PHASE  (GM/CC)  1.799 

THERMAL  CONDUCTIVITY  < J/CM-SEC-DEG  K)  .0016 

THERMAL  DIFFUSIVITY  (CM**2/SE C)  .0006 
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SOLID  PHASE  COMBUSTION  CHARACTERISTICS 


IGNITION  TEMPERATURE  (DEG  K)  300. 

CHEMICAL  ENERGY  (J/GM)  2489. 


MAX  PRESSURE 
(MPA) 


ADD  CONSTANT 
( CM/SEC ) 


PRE-EXPONENT 
( CM/SEC -MPA  *  *BN ) 


EXPONENT 


.  52 
690 . 00 


0  . 
0. 


2. 508 
2 . 007 


.  4620 
.  1330 


GRAIN  GEOMETRY 


EXTERNAL  DIAMETER  (CM)  .3 

LENGTH  (CM)  0. 

DIAMETER  OF  PERFORATIONS  (CM)  0. 

NUMBER  OF  PERFORATIONS  0. 


PROPERTIES  OF  PROJECTILE 


PROJECTILE  MASS  (KG)  46.72 

NUMBER  OF  ENTRIES  IN  BORE  RESISTANCE  TABLE  7 
RESISTANCE  LAW  NUMBER  1 
NUMBER  OF  FILLER  ELEMENTS  0 


BORE  RESISTANCE  DATA 


PROJECTILE  TRAVEL 
(CMS) 

0.000 
1.016 
2 . 540 
3. 937 
5  .  207 
11.430 
520. 700 


RESISTIVE  PRESSURE 
(MPA) 

1 . 72 

23. 10 

34 . 10 
25.00 
22. 40 
17.20 
10. 30 


BAG  FLOW  RESISTANCE  DATA 


TYPE 

INIT  FRICTION 
FACTOR 

RUPTURE  STRESS 

(MPA) 

RUPTURE  INTERVAL 
(MSEC ) 

1 

.01 

.  30 

0 . 0 

2 

101 . 00 

.  60 

0 . 0 

3 

101 .00 

.  60 

2.0 

4 

101 .00 

.  30 

0 . 0 
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DATA  TO  DESCRIBE  REACTIVITY  OF  BAG 
ENERGY  RELEASED  DURING  DECOMPOSITION  (J/GM) 


SUBSTRATE 

2489. 


1 


BAG  SUBSTRATE  DISCHARGE  CHARACTERISTICS 

TIME 

RATE  OF  DISCHARGE 

(MSEC) 

(GM/CM**2-SEC) 

0.0 

2 . 62 

0.1 

26.20 

30.0 

26 . 20 

LOCATION  OF 

POINTS 

FOR  PRESSURE  SUMMARY  TABLE 

AXIAL  LOCATION 
(CMS) 

HALL 

(0 )  OR  AXIS  ( 1 ) 

.  25 

0 

44.45 

0 

88. 27 

0 
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DISTRIBUTION  LIST 


No.  Of 

Copies  Organization 

12  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22314 

1  Office  of  the  Under  Secretary 

of  Defense 

Research  &  Engineering 
ATTN:  R.  Thorkildsen 
Washington,  DC  20301 

1  HQDA/SAUS-OR,  D.  Hardison 

Washington,  DC  20301 

1  HQDA/DAMA-ZA 
Washington,  DC  20310 

2  HQDA,  DAMA-CSM,  A.  German 

E.  Lippi 

Washington,  DC  20310 

1  HQDA/ SARD A 

Washington,  DC  20310 

1  Commander 

US  Army  War  College 
ATTN:  Library-FF229 
Carlisle  Barracks,  PA  17013 

1  Ballistic  Missile  Defense 

Advanced  Technology  Center 
P.  0.  Box  1500 
Huntsville,  AL  35804 

1  Chai rman 

D0D  Explosives  Safety  Board 
Room  856-C 
Hoffman  Bldg.  1 
2461  Eisenhower  Avenue 
Alexandria,  VA  22331 


No.  Of 

Copies  Organization 

3  Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMD-ST 

DCRSF-E,  Safety  Office 
DRCDE-DW 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

14  Commander 

US  Army  Armament  R&D  Command 
ATTN:  DRDAR-TSS  (2  cys) 
DRDAR-TDC 

D.  Gyorog 
DRDAR-LCA 

K.  Russell 

A.  Moss 

J .  Lannon 

A.  Beardell 

D.  Downs 

S.  Einstein 

L.  Schlosberg 
S.  Westley 

S.  Bernstein 
P.  Kemmey 
C .  Heyman 
Dover,  NJ  07801 

9  US  Army  Armament  R&D  Command 

ATTN:  DRDAR-SCA*  L.  Stiefel 

B.  Brodman 

DRDAR-LCB-I,  D.  Spring 
DRDAR-LCE ,  R.  Walker 
DRDAR-LCU-CT 

E.  Barrieres 

R.  Davit t 
DRDAR-LCU-CV 

C. Mandala 
E.  Moore 

DRDAR-LCM-E 

S.  Kaplowitz 
Dover,  NJ  07801 
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No.  Of 
Copies 


DISTRIBUTION  LIST 


No.  Of 

Organization  Copies  Organization 


5  Commander 

US  Army  Armament  R&D  Command 
ATTN:  DRDAR-OAR, 

J .  Rutkowski 

G.  Allen 
J .  Donner 
P.  Serao 
D ,  Adams 

Dover,  NJ  07801 

5  Project  Manager 

Cannon  Artillery  Weapons 
Systems 

ATTN:  DRCPM-CAWS, 

F .  Menke 
DRCPM-CAWS-WS 

H.  Noble 
DRCPM-CAWS- SI 

M.  Fisette 
DRCPM-CAWS- AM 
R.  DeKleine 
H.  Hassman 
Dover,  NJ  07801 

3  Project  Manager 

Munitions  Production  Base 
Modernization  and  Expansion 
ATTN:  DRCPM-PMB,  J.  Ziegler 

M.  Lohr 

A.  Sikloski 
Dover,  NJ  07801 

3  Project  Manager 

Tank  Main  Armament  System 
ATTN:  DRCPM-TMA,  D.  Appling 

DRCPM-TMA- 105 
DRCPM-TMA- 120 
Dover,  NJ  07801 

4  Commander 

US  Army  Armament  R&D  Command 
ATTN:  DRDAR-LCW-A 

M. Salsbury 
DRDAR-LCS 

DRDAR-LCU,  A.  Moss 
DRDAR-LC ,  J.  Frasier 
Dover,  NJ  07801 


5  Commander 

US  Army  Armament  Materiel 
Readiness  Command 
ATTN:  DRDAR-LEP-L,  Tech  Lib 

DRSAR-LC ,  L.  Ambrosini 
DRSAR-ICR,  G.  Cowan 
DRSAR-LEM,  W.  Fortune 
R.  Zastrow 

Rock  Island,  IL  61299 

1  Commander 

US  Army  Watervliet  Arsenal 
ATTN:  SARWV-RD,  R.  Thierry 

Watervliet,  NY  12189 

1  Director 

US  Army  ARRADCOM  Benet 
Weapons  Laboratory 
ATTN:  DRDAR-LCB-TL 

Watervliet,  NY  12189 

1  Commander 

US  Army  Aviation  Research 
and  Development  Command 
ATTN:  DRDAV-E 

4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120 

1  Commander 

US  Army  TSARCOM 
4300  Goodfellow  Blvd 
St.  Louis,  MO  63120 


1  Director 

US  Army  Air  Mobility  Research 
And  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 
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DISTRIBUTION  LIST 


No .  Of 

Copies  Organization 

1  Commander 

US  Army  Communications 
Research  and  Development 
Command 

ATTN:  DRDCO-PPA-SA 

Fort  Monmouth,  NJ  07703 

1  Commander 

US  Army  Electronics  Research 
and  Development  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 

Fort  Monmouth,  NJ  07703 

1  Commander 

US  Army  Harry  Diamond  Lab* 
ATTN:  DELHD-TA-L 
2800  Powder  Mill  Road 
Adelphi ,  MD  20783 

2  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-R 

DRSMI-YDL 

Redstone  Arsenal,  AL 

1  Commander 

US  Army  Natick  Research  and 
Development  Command 
ATTN:  DRXRE ,  D.  Sieling 

Natick,  MA  01762 

1  Commander 

US  Army  Tank  Automotive 

Research  and  Development 
Command 

ATTN:  DRDTA-UL 

Warren,  MI  48090 

1  US  Army  Tank  Automotive 

Research  and  Development  Cmd 
ATTN:  DRSTA-CG 

Warren,  MI  48090 


No.  Of 

Copies  Organization 

1  Project  Manager 
Improved  TOW  Vehicle 
ATTN:-  DRCPM-ITV 

US  Army  Tank  Automotive 
Research  &  Development 
Command 

Warren,  MI  48090 

2  Program  Manager 

Ml  Abrams  Tank  System 
ATTN :  DRCPM-GMC-S A , 

J.  Roossien 
Warren,  MI  48090 

1  Project  Manager 

Fighting  Vehicle  Systems 
ATTN:  DRCPM-FVS 

Warren,  MI  48090 

1  Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL,  Tech  Lib 
White  Sands  Missile  Range, 

NM  88002 

1  Project  Manager 

M-60  Tank  Development 
ATTN:  DRCPM-M60TD 

Warren,  MI  48090 

1  Commander 

US  Army  Training  &  Doctrine 
Command 

ATTN:  ATCD-MA/  MAJ  Williams 

Fort  Monroe,  VA  23351 

2  Commander 

US  Army  Materials  and 
Mechanics 
Research  Center 
ATTN :  DRXMR-ATL 

Tech  Library 
Watertown,  MA  02172 
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DISTRIBUTION  LIST 


No.  Of 

Copies  Organization 

1  Commander 

US  Army  Research  Office 
ATTN :  Tech  Library 
P.  0.  Box  12211 
Research  Triangle  Park,  NC 
27709 

1  Commander 

US  Army  Mobility  Equipment 
Research  &  Development 
Command 

ATTN:  DRDME-WC 
Fort  Belvoir,  VA  22060 

1  Commander 

US  Army  Logistics  Mgmt  Ctr 
Defense  Logistics  Studies 
Fort  Lee,  VA  23801 

2  Commandant 

US  Army  Infantry  School 
ATTN:  Infantry  Agency 

Fort  Benning,  GA  31905 

1  US  Army  Armor  &  Engineer 

Board 

ATTN:  STEBB-AD-S 
Fort  Knox,  KY  40121 

1  Commandant 

US  Army  Aviation  School 
ATTN:  Aviation  Agency 
Fort  Rucker,  AL  36360 

1  Commandant 

Command  and  General  Staff 
College 

Fort  Leavenworth,  KS  66027 

1  Commandant 

US  Army  Special  Warfare 
School 

ATTN:  Rev  &  Tng  Lit  Div 
Fort  Bragg,  NC  28307 

1  Commandant 

US  Army  Engineer  School 

ATTN :  ATSE-CD 

Ft.  Belvoir,  VA  22060 


No.  Of 

Copies  Organization 

1  Commander 

US  Army  Foreign  Science  & 

Te chn o 1 o gy  Center 
ATTN:  DRXST-MC-3 
220  Seventh  Street,  NE 
Charlottesville,  VA  22901 

1  President 

US  Army  Artillery  Board 
Ft.  Sill,  OK  73504 

2  Commandant 

US  Army  Field  Artillery 
School 

ATTN:  ATSF-C0-MW,  B.  Willis 
Ft.  Sill,  OK  73503 

3  Commandant 

US  Army  Armor  School 
ATTN:  ATZK-CD-MS/M. 

Falkovitch  Armor  Agency 
Fort  Knox,  KY  40120 

1  Chief  of  Naval  Materiel 

Department  of  the  Navy 
ATTN:  J.  Amlie 
Washington,  DC  20360 

1  Chief  Naval  Research 

ATTN:  Code  473,  R.  S.  Miller 

800  N.  Ouincy  Street 
Arlington,  VA  22217 

2  Commander 

US  Naval  Sea  Systems  Command 
ATTN:  SEA-62R2,  J.  W.  Murrin 

R.  Beauregard 
Washington,  DC  20362 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  NAIR-954-Tech  Lib 
Washington,  DC  20360 
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1  Strategic  Systems  Project 

Office 

Dept,  of  the  Navy 
Room  901 

ATTN:  J.  F.  Kincaid 

Washington,  DC  20376 

1  Assistant  Secretary  of  the 

Navy  (R,  E,  and  S) 

ATTN:  R.  Reichenbach 

Room  5E787 
Pentagon  Bldg. 

Washington,  DC  20350 

1  Commander 

Naval  Research  Lab 

Tech  Library 
Washington,  DC  20375 

5  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  G33,  J.  L.  East 

D.  McClure 
W.  Burrell 
J.  Johndrow 
Code  DX-21  Tech  Lib 
Dahlgren,  VA  22448 

2  Commander 

US  Naval  Surface  Weapons 
Center 

ATTN:  J.  P.  Consaga 

C.  Gotzmer 

Indian  Head,  MD  20640 

4  Commander 

Naval  Surface  Weapons  Center 
ATTN:  S.  Jacobs /Code  240 

Code  730 

K.  Kim/Code  R-13 
R.  Bernecker 
Silver  Spring,  MD  20910 

2  Commanding  Officer 

Naval  Underwater  Systems 
Center 

Energy  Conversion  Dept. 

ATTN:  CODE  5B331,  R.  S.  Lazar 
Tech  Lib 

Newport,  RT  02840 


4  Commander 

Naval  Weapons  Center 
ATTN:  Code  388,  R.  L.  Derr 

C.  F.  Price 
T.  Boggs 
Info.  Sci .  Div. 

China  Lake,  CA  93555 

2  Superintendent 

Naval  Postgraduate  School 
Dept,  of  Mechanical 
Engineering 
ATTN:  A.  E.  Fuhs 

Code  1424  Library 
Monterey,  CA  93940 

6  Commander 

Naval  Ordnance  Station 
ATTN:  P.  L.  Stang 

C •  Smi th 
S.  Mitchell 

C.  Christensen 

D.  Brooks 
Tech  Library 

Indian  Head,  MD  20640 

1  HQ  AFSC 

Andrews  AFB 
Washington,  DC  20331 

1  Program  Manager 

AFOSR 

Directorate  of  Aerospace 
Sciences 

ATTN:  L.  H.  Caveny 

Bolling  AFB,  DC  20332 

6  AFRPL  (DYSC) 

ATTN:  D.  George 

J.  N.  Levine 
B.  Goshgarian 
D.  Thrasher 
N.  Vander  Hyde 
Tech  Library 
Edwards  AFB,  CA  93523 
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AFFTC 

ATTN:  SSD-Tech  Lib 

Edwards  AFB,  CA  93523 

AFATL 

ATTN:  DLYV 
Eglin  AFB,  FL  32542 

AFATL/DLDL 

ATTN:  0.  K.  Heiney 

Eglin  AFB,  FL  32542 


1  AVCO -Everett  Rsch  Lab  Div 

ATTN:  D.  Stickler 

2385  Revere  Beach  Parkway 
Everett,  MA  02149 

2  Calspan  Corporation 
ATTN:  E.  B.  Fisher 

Tech  Library 
P.  0.  Box  400 
Buffalo,  NY  14225 


ADTC 

ATTN:  DLODL  Tech  Lib 

Eglin  AFB,  FL  32542 

AFFDL 

ATTN:  TST-Lib 

Wright-Patterson  AFB,  OH 
45433 


1 


1 


Foster  Miller  Associates 
ATTN:  A.  Erickson 

135  Second  Avenue 
Waltham,  MD  02154 


Hercules,  Inc. 
Bacchus  Works 
ATTN:  K.  P.  McCarty 

Magna,  UT  84044 


HQ  NASA 

600  Independence  Avenue,  SW  i 

ATTN:  Code  JM6,  Tech  Lib. 

Washington,  DC  20546 

NASA/Lyndon  B.  Johnson  Space 
Center 

ATTN:  NHS-22 ,  Library  i 

Section 

Houston,  TX  77058 

Aerodyne  Research,  Inc. 

Bedford  Research  Park 
ATTN:  V.  Yousef ian 

Bedford,  MA  01730  1 

Aerojet  Solid  Propulsion  Co. 

ATTN:  P.  Micheli 

Sacramento,  CA  95813 

Atlantic  Research  Corporation 
ATTN:  M.  K.  King 
5390  Cheorokee  Avenue 
Alexandria,  VA  22314 


General  Applied  Sciences  Lab,  Inc. 
ATTN:  J.  Erdos 

Merrick  &  Stewart  Avenues 
Westbury  Long  Island,  NY 
11590 

General  Electric  Company 
Armament  Systems  Dept. 

ATTN:  M.  J.  Bulman, 

Room  1311 
Lakeside  Avenue 
Burlington,  VT  05412 

Hercules  Powder  Co. 

Allegheny  Ballistics 
Laboratory 
ATTN:  R.  B.  Miller 

P.  0.  Box  210 
Cumberland,  MD  21501 
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1  Hercules,  Inc. 

Eglin  Operations  2 

AFATL  DLDL 

ATTN:  R.  L.  Simmons 

Eglin  AFB ,  FL  32542 

1  IITRI 

ATTN:  M.  J.  Klein 

10  W.  35th  Street  1 

Chicago,  IL  60616 

2  Lawrence  Livermore  Laboratory 

ATTN:  M.  S.  L-355, 

A.  Buckingham  1 

M.  Finger 
P.  0.  Box  808 
Livermore,  CA  94550 


1  Olin  Corporation 

Badger  Army  Ammunition  Plant  i 

ATTN:  R.  J.  Thiede 

Baraboo ,  WI  53913 


1  Olin  Corporation 

Smokeless  Powder  Operations  3 

ATTN:  R.  L.  Cook 
P.  0.  Box  222 
ST.  Marks,  FL  32355 


1  Paul  Gough  Associates,  Inc. 

ATTN:  P.  S.  Gough 

P.  0.  Box  1614  2 

Portsmouth,  NH  03801 

1  Physics  International  Company 

2700  Merced  Street 
Leandro,  CA  94577 


Organization 


Rockwell  International 
Rocket dyne  Division 
ATTN:  BA08  J.  E.  Flanagan 

J.  Grey 

6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

Science  Applications,  INC. 
ATTN:  R.  B.  Edelman 

23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

Scientific  Research  Assoc., 
Inc . 

ATTN:  H.  McDonald 

P.  0.  Box  498 
Glastonbury,  CT  06033 

Shock  Hydrodynamics,  Inc. 
ATTN:  W.  H.  Andersen 

4710-16  Vineland  Avenue 
North  Hollywood,  CA  91602 

Thi oko 1  Co  rp  o  r a  t i on 
Huntsville  Division 
ATTN:  D.  Flanigan 

R.  Glick 
Tech  Library 
Huntsville,  AL  35807 

Thiokol  Corporation 
Wasatch  Division 
ATTN:  J.  Peterson 

Tech  Library 
P.  0.  Box  524 
Brigham  City,  UT  84302 


2  Thiokol  Corporation 

Elkton  Division 
ATTN:  R.  Biddle 

Tech  Lib. 
Elkton,  MD  21921 

1  Pulsepower  Systems,  Inc. 

ATTN:  L.  C.  Elmore 

815  American  Street 
San  Carlos,  CA  94070 


1  Princeton  Combustion  Research 

Lab .  ,  Inc . 

ATTN:  M.  Summerfield 

1041  US  Highway  One  North 
Princeton,  NJ  08540 
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2  United  Technologies 

Chemical  Systems  Division 
ATTN:  R.  Brown 

Tech  Library 
P.  0.  Box  358 
Sunnyvale,  CA  94086 

1  Universal  Propulsion  Company, Inc 

ATTN:  H.  J.  McSpadden 

Black  Canyon  Stage  1,  Box  1140 
Phoenix,  AZ  85029 

1  Southwest  Research  Institute 

ATTN:  W.  H.  McLain 

8500  Culebra  Road 
San  Antonio,  TX  98228 


1  Battelle  Memorial  Institute 

ATTN:  Tech  Library 

505  King  avenue 
Columbus,  OH  43201 

1  Brigham  Young  University 

Dept,  of  Chemical  Engineering 
ATTN:  M.  Becks tead 
Provo,  UT  84601 


1  California  Institute  of  Tech 

204  Karman  Lab 

Main  Stop  301-46 
ATTN:  F.  E.  C.  Culick 

1201  E.  California  Street 
Pasadena,  CA  91125 

1  California  Institute  of  Tech 

Jet  Propulsion  Laboratory 
ATTN:  L.  D.  Strand 
4800  Oak  Grove  Drive 
Pasadena,  CA  91103 

1  University  of  Illinois 

Dept  of  Mech  Eng 
ATTN:  H.  Krier 

144  MEB,  1206  W.  Green  St. 
Urbana,  IL  61801 


1  University  of  Massachusetts 

Dept,  of  Mechanical 
Engineering 
ATTN:  K.  Jakus 

Amherst,  MA  01002 

1  University  of  Minnesota 

Dept,  of  Mechanical 
Engineering 
ATTN:  E.  Fletcher 

Minneapolis,  MN  55455 

1  Case  Western  Reserve 

Uni versi ty 

Division  of  Aerospace 
Sciences 
ATTN:  J.  Tien 
Cleveland,  OH  44135 

3  Georgia  Institute  of  Tech 

School  of  Aerospace  Eng. 

ATTN:  B.  T.  Zinn 

E.  Price 
W.  C.  Strahle 
Atlanta,  GA  30332 

1  Institute  of  Gas  Technology 

ATTN:  D.  Gidaspow 

3424  S.  State  Street 
Chicago,  IL  60616 

1  Director 

Applied  Physics  Lab/  Chemical 
Propulsion  Info  Agency 
The  Johns  Hopkins  University 
ATTN:  T.  Christian 

Johns  Hopkins  Road 
Laurel,  MD  20707 

1  Massachusetts  Institute  of 

Tech 

Dept  of  Mechanical 
Engineeri ng 
ATTN:  T.  Toong 

Cambridge,  MA  02139 
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1  Pennsylvania  State  University 

Applied  Research  Lab 
ATTN:  G.  M.  Faeth 

P.  0.  Box  30 
State  College,  PA  16801 

1  Pennsylvania  State  University 

Dept.  Of  Mechanical 
Engineering 
ATTN:  K.  Kuo 

University  Park,  PA  16802 

1  Purdue  University 

School  of  Mechanical 
Engineering 
ATTN:  J.  R.  Osborn 

TSPC  Chaffee  Hall 
West  Lafayette,  IN  47906 

1  Rensselaer  Polytechnic  Inst. 

Department  of  Mathematics 
Troy,  NY  12181 

1  Rutgers  University 

Dept,  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  S.  Temkin 

University  Heights  Campus 
New  Brunswick,  NJ  08903 

1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Tech  Library 

333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

1  Stevens  Institute  of 

Technology 
Davidson  Laboratory 
ATTN:  R.  McAlevy ,  III 

Hoboken,  NJ  07030 

2  Los  Alamos  National  Lab 
ATTN:  T.  D.  Butler,  MS  B216 

M.  Division,  B.  Craig 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 


1  University  of  Southern 

California 

Mechanical  Engineering  Dept. 
ATTN:  0HE200 ,  M.  Gerstein 

Los  Angeles ,  CA  90007 

2  University  of  Utah 

Dept,  of  Chemical  Engineering 
ATTN:  A.  Baer 

G.  Flandro 

Salt  Lake  City,  UT  84112 

1  Washington  State  University 

Dept,  of  Mechanical 
Engineering 
ATTN:  C.  T.  Crowe 

Pullman,  WA  99164 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  H.  Cohen 
Cdr ,  USATECOM 

ATTN:  DRSTE-TO-F 

STEAP-MT ,  S.  Walton 
G.  Rice 
D.  Lacey 
C.  Herud 

Dir,  HEL 

ATTN:  J.  Weisz 

Dir,  USACSL ,  Bldg.  E3516,  EA 
ATTN:  DRDAR-CLB-PA 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports . 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name  : _ _ _ 

Telephone  Number: _ _ _ _ _ 

Organization  Address:  _ _ _ _ _ 


